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Tomonaga-Luttinger liquid correlations and Fabry-Perot interference in conductance 
and finite-frequency shot noise in a single-walled carbon nanotube 
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We present a detailed theoretical investigation of transport through a single-walled carbon nan- 
otube (SWNT) in good contact to metal leads where weak backscattering at the interfaces between 
SWNT and source and drain reservoirs gives rise to electronic Fabry-Perot (FP) oscillations in 
conductance and shot noise. We include the electron-electron interaction and the finite length of 
the SWNT within the inhomogeneous Tomonaga-Luttinger liquid (TLL) model and treat the non- 
equilibrium effects due to an applied bias voltage within the Keldysh approach. In low-frequency 
transport properties, the TLL effect is apparent mainly via power-law characteristics as a function 
of bias voltage or temperature at energy scales above the finite level spacing of the SWNT. The 
FP-frequency is dominated by the non-interacting spin mode velocity due to two degenerate sub- 
bands rather than the interacting charge velocity. At higher frequencies, the excess noise is shown 
to be capable of resolving the splintering of the transported electrons arising from the mismatch of 
the TLL-parameter at the interface between metal reservoirs and SWNT. This dynamics leads to 
a periodic shot noise suppression as a function of frequency and with a period that is determined 
solely by the charge velocity. At large bias voltages, these oscillations are dominant over the ordi- 
nary FP-oscillations caused by two weak backscatterers. This makes shot noise an invaluable tool 
to distinguish the two mode velocities in the SWNT. 

PACS numbers: 73.63.-b, 71.10.Pm, 72.10.-d, 72.70.-fm, 73.63.Fg 
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I. INTRODUCTION 

The study of one-dimensional (ID) electronic systems 
has attracted much interest due to its unique propertiesi. 
In ID, electron-electron (e-e) interaction cannot be ne- 
glected anymore but changes the physical properties 
drastically unlike in higher dimensional metals which are 
described successfully by the Fermi liquid theory. More 
specifically, the notion of quasiparticle excitations com- 
pletely breaks down in ID, and the low-energy excitations 
are collective charge and spin modes travelling at differ- 
ent speeds, a phenomenon known as spin-charge separa- 
tion. 

The low-energy properties of ID metals have been 
investigated successfully within the framework of the 
Tomonaga-Luttinger liquid (TLL) theory^'^. Recently, 
a renewed interest in ID systems has emerged due to 
the possibility to fabricate ideal ID conductors like car- 
bon nanotubes or semiconductor quantum wires. In- 
deed, characteristic predictions of the TLL-model like 
power-law renormalized conductanc o'^i^i^ or spin-charge 
separation— have been confirmed in the tunneling regime 
where the ID system is well separated from the higher- 
dimensional reservoirs. Only recently, transport experi- 
ments through single-walled carbon nanotubes (SWNTs) 
with an average conductance close to the theoretical max- 
imum of Go — Ae^/h, h is the Planck constant and e is the 
electron charge, have been achievecl2i2ii£. On the theoret- 
ical side, Pega et al. have calculated the zero temperature 
conductance for the model of a SWNT in good contact to 
two metal reservoirs and found that the Fabry-Perot (FP) 
type of interference due to phase coherent motion within 



the SWNT is modified by electron-electron interactionii. 
To our knowledge, current noise has not yet been calcu- 
lated in this regime of weak backscattering including the 
FP-interference between two barriers, see Fig. I, as well 
as two spinfull bands which seems crucial to understand 
existing shot noise experiments in SWNT— where some 
weak backscattering at the SWNT-metal reservoir inter- 
face cannot be avoided. 

Shot noise is sensitive to temporal correlations of the 
current and thus provides additional information about 
dynamical processes inside the conductor not accessible 
in conductancei^. In particular, noise is sensitive to the 
elementary excitations of the system. In the edge states 
of the fractional quantum Hall effect regime, a chiral TLL 
is realized, where right- and left-going particles are lo- 
cated at different edges of the sample. The fractional 
charge ge, with g the TLL-parameter, has been mea- 
sured in low-frequency shot noise^^ in agreement with 
theoryii. Shot noise measurements in SWNT are very 
recen1i^&^^i^ and no quantitative analysis of shot noise 
measurements in the TLL-regime have been reported so 
far. In a SWNT right- and left-moving electrons coexist 
in the same channel, and consequently electrons can scat- 
ter at the interface between the TLL-system and the non- 
interacting reservoirs. Therefore, the physics is expected 
to be quite different from its chiral counterpart. One pos- 
sibility to model the finite size effect and the influence of 
the reservoirs is to use the inhomogeneous TLL-model 
where the interaction parameter changes from 17 = I in 
the reservoirs to g < I in the interacting regioni^iiSiiS. 
Within this model it has been found that the fractional 
charge ge of the TLL cannot be simply extracted from 



the ratio between shot noise and backscattered current. 
It is rather the stable charge e of the reservoir carriers to 
which shot noise is sensitive at low frequencies. This has 
been concluded for a single-channel TLL with spin sub- 
jected to a random backscattering potential^" and for a 
single channel spinless TLL with a single impurity within 
the wireSi*2SiS. We reach here the same conclusion in 
the specific FP-setup of Fig. 1. Despite the lack of a di- 
rect measurement of the fractional charge through low- 
frequency noise properties, the low- frequency shot noise 
is sensitive to interaction since the backscattering off the 
barriers is energy dependent leading to power-law depen- 
dent noise S and Fano factor F — S/el, where / denotes 
the average current. 

Recently, it became possible to measure also high- 
frequency nois o^^i^^ . This opens up a way to explore in- 
teraction related effects in an extended parameter range. 
As shown in Refs. 22,23,26, the high-frequency noise be- 
comes sensible to the momentum- conserving reflections 
of charge excitations due to the mismatch of g at the in- 
terface between the SWNT and metal reservoirs which 
allows to extract further information about g not con- 
tained in low-frequency transport properties. These mul- 
tiple reflections are even present without any physical 
scatterer— , but are only resolved in transport for fre- 
quencies on the order of the interacting level spacing, i.e. 
huj ^ hvF/2Lg where Vp is the Fermi velocity and L is 
the length of the interacting region. However, the situa- 
tion is different once an impurity is included in the sys- 
tem. Electron waves can be scattered at the impurity site 
and interfere with the transmitted part which is partially 
backscattered at the interface due to the inhomogene- 
ity of g which leads also to oscillations with frequency 
v^/^Lg as a function of bias voltage. This point was 
noted in several worksiit^iSLSS. In the experimentally 
relevant case of a SWNT with two impurities, this inter- 
action induced interference is masked by the usual FP- 
oscillations due to two scatterers naturally formed at the 
interface between the SWNT and metal reservoirs. Since 
the SWNT has three non-interacting modes due to spin 
and subband degeneracy and only one interacting mode 
of the total charge carrying the information about g, any 
oscillation in the bias voltage dependence of conductance 
or noise is dominated by the non-interacting spin mode 
frequency Up/i. However, as pointed out in Ref. 11, ap- 
plying a gate voltage can decrease the amplitude of the 
ordinary FP-interference. In that case, small oscillations 
with frequency v^/iLg remain. The TLL-parameter g 
is presumably only weakly dependent on gate voltage^-^. 
However, in general, applying a gate voltage can influence 
g in a TLL due to screening by the gate electrode^*^^. We 
find now that noise as a function of frequency to and bias 
voltage V is capable of clearly discriminating the two os- 
cillation periods of collective modes present in the SWNT 
without changing the gate voltage. At high bias voltages 
(eV ^ huj^ hvp/L), we find that the frequency dependent 
excess noise shows oscillations dominated by the charge- 
mode frequency v-p/2Lg^ whereas the bias voltage de- 



pendence exhibits the FP-oscillations dominated by the 
non-interacting spin- mode frequency v^/ L whose ampli- 
tude is modulated by lu. This clearly distinguishes the 
charge plasmon resonance induced by the finite length 
L of the interacting region (SWNT) from the more con- 
ventional FP-interference due to two barriers. The finite 
frequency noise therefore could be used to extract both 
frequency scales which allows us to extract g without the 
knowledge of any system parameters like the position of 
an impuritySSiS^ or the fitting to a power-law^. This 
is highly anticipated since power-laws can also originate 
from environmental effects (dynamical Coulomb block- 
ade) in the same functional form^. 

The organization of the paper is as follows: in Section 
II we introduce the TLL-model of a SWNT with spatially 
inhomogeneous TLL-parameter taking into account the 
effects of the non-interacting source and drain electrodes. 
We then discuss the inclusion of two weak backscattering 
potentials situated at the interfaces between metal elec- 
trodes and SWNT. In Section III we introduce the gen- 
eral framework of a Keldysh functional integral approach 
to treat the non-equilibrium effects due to an applied bias 
voltage. In Section IV we present the dc conductance to 
leading order in the backscattering thereby extending the 
result of Ref. 11 to finite temperatures^i. We also give 
some asymptotic analytical results showing the relevant 
power-laws and provide numerical results for the generic 
case. Section V is devoted to the current noise where 
we discuss the low-frequency noise, Fano factor and the 
general frequency dependence. The details of the calcu- 
lations are presented in the appendices. Sections IV and 
V close with a discussion of the physical interpretation 
of the results. We set ?i = 1 in intermediate steps but 
restore fi in final results. 



II. MODEL FOR SWNT COUPLED TO METAL 
RESERVOIRS 

We consider electrons in a SWNT subjected to a re- 
pulsive Coulomb interaction potential parametrized by 
A > with Hamiltonian density'^^i^'^ 



Hs 






(1) 



where ptot{x) = X^Li Es=T,i('^LsV'ffis+V'LsV'Ljs) is the 
total charge density and i = 1,2 denotes the two bands 
that cross the Fermi level. The slow- varying parts of the 
field-operators for left(L) and right(R) moving electrons 
can be expressed in terms of bosonic fields as 



satisfy 



^/2^ 



(2) 



which 

[cj)is{x),9js'{x')\ 



the commutation relation 

i{'K /2)6ij5ss'S.gTi{x — x'). This 



relation implies that Ilis{x) — ~{l/TT)dx4>is{x) is the 
conjugate momentum to 9is{x). In Eq. (|2J) we have 
introduced a short-distance cut-ofF A which is on the 
order of the lattice spacing^''. To proceed, it is useful to 
define new fields for total charge(spin) and charge(spin) 
imbalance between the two bands. We define charge(c) 
and spin((T) bosonic fields via 9ic = {9i-\ + 9ii)/^/2 and 
(^ia = (^JT ~ ^ii)/V2 and further the symmetric (+) and 
antisymmetric(-) combinations 0±^ — (^i^ ± 6'2^)/\/2, 
^ = c,(T and similar for (/i-fields. We obtain four labek^^: 
a — {I = +p,2 = +a,3 = — p, 4 = — c}. In this 
new basis the Hamiltonian H^wnt — J dxHswNT for the 
SWNT incorporating the reservoirs becomes 



H, 



2tt 



dx 



{dx4>if 



2-z:x(^-^i^' 



g'^{x) 



r 

+ ?E / ^^ [{d.<^a? + {dJaf] ■ (3) 

^ a=2 •' 

The velocity of the collective charge excitations in the 
SWNT is Vc — Vp/g which is renormalized due to re- 
pulsive e-e interaction in the nanotube. Since the in- 
teraction potential strength A couples only to the total 
charge density, only the charge sector a = 1 is modified 
by the TLL-parameter g = [vp/{v^ + {8\/tt))]^/^. We 
assume g{x) = g < 1 in the SWNT and g{x) = 1 in the 
reservoirs. The inhomogeneity of g reflects the finite size 
of the nanotube. The abrupt change of g at the inter- 
faces between metal reservoirs and SWNT is considered 
to be a good approximation to a smooth transition of g 
as long as the real length over which g changes is much 
smaller than the typical wavelengths of the excitations in 
the TLL, but larger than the Fermi wavelength or lattice 
spacing22,. The relation of the bosonic fields in Eq. Q to 
physical quantities can be examined by looking at prod- 
ucts of fermion operators. Using the relation for normal 
ordered densities ^^/^(a;) =: ip]^^^^^{x)ipii/L,,{x) := 
±dx{(t)is{x) ± 9is{x))/2Tr, we obtain, e.g. for the total 
charge density, ptot{x) = {2/T:)dx9i{x). Of particular in- 
terest is the operator for the charge current. From the 
continuity equation we obtain I{x,t) = —e{2/T:)9i{x,t). 
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FIG. 1: The Fabry-Perot double barrier device: The backscat- 
tering with bare amplitude strengths ui and U2 of electrons 
primarily takes part at the SWNT-metal reservoir interfaces 
where the TLL parameter g changes from g=l in the leads to 
g < 1 in the nanotube. 

The backscattering off impurities is assumed to be 
weak and mainly happening at the two metal contact- 



SWNT interfaces which separate the nanotube from the 
reservoirs. The form of the backscattering Hamiltonian 
is given as 



^- = EE E 

m=l i,j=l s=T,i 



-j^^H-D^+^A. 



^'4'\is{Xm)ll^R]s{Xm) + h-C. 
2 

= X! E "™ exp{i[6'i,„ + s92r. 

m,^,^ — 1 S— ±1 
-H(-1)'+My(03m+S04™) 

+(-l)"+iA,,]}+h.c. 



(4) 



In Eq. Q) we have used 9am = 9a{xm) and similarly for 
't'am with xi_2 = T-^/2 denoting the positions of the two 
barriers. We have further defined u^ = u^/27rA which 
are real valued and have the dimension of energy. The 
backscattering phase for the scattering of a right-moving 
electron with band-index j to a left-moving electron with 
band-index i is denoted by (— l)™+^Aij. Its dependence 
on the contact label m — 1,2 refiects the mirror sym- 
metry of the two SWNT-metal reservoir interfaces with 
respect to a; = 0. We further assume that these phases 
are energy independent. 

Next, we discuss the inclusion of a gate voltage Vg 
which gives rise to a Hamiltonian density proportional to 
the total charge density H' oc ptot V^ = (2/7r)(9a;6'i)V^. 
This linear term in the Hamiltonian can be eliminated by 
performing the linear shift^^ 9i ^ 9i ~ VgX which leaves 
the quadratic Hamiltonian Eq. Q unchanged (up to an 
irrelevant constant) but changes i?bs where 9i is replaced 
by 6*1 — VgX. Note that applying a gate voltage induces a 
shift of the Fermi level in the SWNT and metal contacts. 



III. THE TRANSPORT THEORY 

In this section we derive the general framework for cal- 
culating the current and current noise in non-equilibrium 
within the Keldysh functional approach. We start with 
the system Hamiltonian H — i^swNx + ^bs and treat _ffbs 
as the perturbation. The average of an observable O is 
{d{t)) = Ti[pd{t)] where d{t) = e»^(*-*o)Oe-^-f^(*-*«), p 
is the density matrix at time to before -ffbs is switched on 
and Tr means trace. The non-equilibrium effect caused 
by the bias voltage V can be included in the density ma- 
trix. We assume that before the backscattering Hamil- 
tonian i?bs is turned on (at to — *■ — oo) the system has 
a well defined non-equilibrium state determined by sep- 
arate chemical potentials for left and right-movers kept 
fixed by the chemical potentials of the right and left elec- 
tron reservoirs, respectively. The initial density matrix 



therefore takes on the forr 



Pv 



1 



-PHv 



(5) 



with Hv = ffswNT - ^-rNr - ij-lNl and (3 = l/fceT 
with T the temperature, k^ the Boltzmann constant and 
Zy = Tr[exp{—/3Hv)] the partition function. The equi- 
librium chemical potential is defined as zero (a non-zero 
chemical potential can be taken into account by the gate 
voltage) and A^_r/l — J dxnji/L{x). The bias voltage is 
then related to the chemical potentials of left and right- 
movers via ^jlr/l = ±ey/2. As outlined in Ref. 11, it 
is convenient to apply a unitary transformation Uv such 
that UyHyUv = ^swNT + const.. This transforms the 
bias voltage from the density matrix {pv -^ Po) into the 
backscattering Hamiltonian i^bs which receives a time de- 
pendent phase factor in the interaction picture governed 
by the shift 9i ^ di — eVt. In addition, the unitary 
transformation transforms the observable according to 
O -^ UyOUv- In the case of the current operator this 
leads to the shifliii / -> /q + /. The average current 
can then be written as {I{x,t)) = Iq + {Iv{x,t))o- Here, 
Iq = Ae'^V/h is the ideal current without backscattering 
and Iv{x,t) gives rise to the backscattered current 



lB{x,t) = (/y(x,t))o 

- - -iT.rf dt'H^^(,t')\ 

= {TKlK{x,t)e - 



(6) 



Here, we have introduced the Keldysh current operator 
lK{x,t) = (1/2) X]r-^'^('^'^)' *^^ *™^ ordering operator 
Tx along the Keldysh contour depicted in Fig. 2 and 
r = ± which refers to fields defined on the zb-branch of 
that contour. In Eq. © the time dependence of all op- 
erators is due to ifswNT only and (...)o = Tr[/5o...]. A 
similar procedure can be performed for the noise spec- 
tral density S{x,uj) = J dte''^'^^-*'^ S{x;t,t'), where the 
symmetrized current-current correlator is S{x] t, t') = 
{l/2){{6iix,t),SI{x,t')}) where {...} denotes the anti- 
commutator, 51{x,t) — I{x,t) — (!) and (...) = Tr[/5y...] 
with the initial density matrix py discussed before. Us- 
ing again formally I(t) = /q + Iv{t) we obtain S{x, uj) ~ 

(l/2)/dte*"(*-*')({/y(x,t),/y(x,t')})o - 2TTSii0)ll. To 
lowest order in the backscattering, the (5-function contri- 
bution at zero frequency can be neglected and we can 
therefore write the current-current correlator as 



The time ordered correlation functions can be conve- 
niently calculated by means of a functional integral ap- 
proach discussed next. 



A. The generating functional 

The statistical averages in Eqs. (0 and O are con- 
veniently evaluated in terms of the following generating 
functional Z^ 



a •' 

zSo -iJdt'iH+it') - iIbl(i')]L±^,±. 



X exp 



■eVt' 



-i / dt' / dx'Tjix',t')ei{x',t') 



(8) 



We have performed the rotation to new fields 9^ = 9ai: 
i9a/2 which allows the simple representation Ix{x,t) = 
—e{2/TT)9i{x,t). In Eq. we have introduced a source 
field ri{x, t) which does not have a direct physical meaning 
but is rather a convenient way to produce correlation 
functions via of functional derivatives. The action 5*0 
describes the dynamics induced by Hswnt only and is 
a quadratic form of the phase fields 9a{x, t), 9a{x, t) and 
4>a{x, t), 4>a{x, t). The explicit form of Sq is presented in 
Appendix A. Here, we only give the relevant correlation 
functions 



Cl\x,x'-t)={fK9a{x,t)ea{x',Q))o 

= ]^{{9a{x,t),9a[x',G)})„, 



(9) 



and the retarded functions 



i?f (X, x'-t) ={TK9aix, t)9aix', 0))o 

= -iQ{t){[9a{x,t),9aix',0)])o, (10) 



S{x;t,t') = 

{TKlK{x,t)lK{x,t')e ^ ^Ur^n-— ")o. 

(7) 



and similar for ^Q-correlations. Other combinations like 
{9a{x, t)9a{x', 0))o = (Mx, t)Mx', 0))o - 0. 



branch 



B. Shifted action 



Time 



■ branch 



FIG. 2: The Keldysh contour: Operators are ordered along 
the contour with operators evaluated at later times acting on 
the left of operators evaluated at earlier times. Times on the 
(-I-) branch are always earlier than times on the (-) branch. 



It is advantageous to transform away the linear ?7-term 
in the generating functional Z^ by shifting the 6'i-fields 
such that in the new variables the linear term in 9i gets 
cancelled, whereas 5*0 remains unchanged. Since we have 
to perform such a transformation on the whole action, 
including the backscattering contribution, the ry-source 
field will appear in the backscattering Hamiltonian in- 



stead. This transformation we find to be 



and 



-_L / dx' I dujuje-"^'ri{x',Lo)Cl\x,x']Uj), 



i{x,t) -> Oi{x,t) 



-2^ I '^' 



duj uo e-'"'* 77(2:', uj) i?f (x', x; -w). 



(11) 

Since the action Sq couples (pi and 61 (see Appendix A), 
01 gets also transformed. However, its transformation 
is not needed here since (/)i-terms are absent in Hhs [see 
Eq. Q] which states that the total charge is conserved 
in the backscattering process. In the new variables the 
generating functional becomes 

yr) ^ — ^ / du)u) j dx' j dx" rj{x' .u))* C ^ {x' ,x" ■,Lj)r](x" ,u) 



-^I dt'iHl-H-Jx 



(12) 



where ri{x,uj)* — ri{x,~uj) and we used the abbrevia- 
tion (...)o - na/^[^^'/'a]--exp[iS'o]. The arrow -> 
in Eq. l(T^ depicts the shift of Of via Eq. ((TT|l and 
the effect of the apphed voltages, exphcitly Of^it) -^ 
^LW + (l/27r)/dw/dx'we-'"*r/(a;',w)[Cf^(a:„,x';a;)± 
^Ri^{x',Xm; — w)] — eVt ~ VgXm- The generating func- 
tional Eq. (|12ll is the starting point for calculating any 
order of current-current correlation functions for a gen- 
eral measurement position x. 



IV. DC CURRENT 



C/™ = (STT^/a^je-^^^") y <u^^'cos(l/gL-t-2Ay) 



We have introduced the dimensionless time t = t/t-^ 
with ip = L/vf the non-interacting traversal time of the 
SWNT as well as the dimensionless voltage v = eVU/h. 
Note that the dc current is independent of the measure- 
ment point x and time t. Each backscattering event 
involves a combination of the total charge mode {9i9i- 
correlations) and the three non-interacting modes {OaOa- 
or (/)Q(/)a-correlations, a = 2,3,4). Therefore, Cmm'{t) = 
^mm' ('•) + 3 C^„/ (t) and a similar definition holds for 
R-mm'(*) (C -^ R)- The superscripts I and F refer to 
interacting {g < 1) and free (non-interacting, i.e. g = 1), 
respectively. Eq. H14(l is consistent with the conductance 
formula derived in Ref. 11 up to the scattering phases 
Aij which have been neglected previously. Physically, 
the term in Eq. H14() proportional to C/'" describes the 
incoherent addition of two barriers whereas the term pro- 
portional to C/'^° describes the quantum mechanical inter- 
ference between backscattering events of different bar- 
riers (1 or 2). Note that the interference term can be 
modulated by the gate voltage Vg. In addition, different 
scattering phases Ay for intraband (i — j) and interband 
{i y^ j) processes lead to a reduction of the FP-amplitude 
(see also Ref. 8). A general analytical form of the con- 
ductance seems difficult to derive and we have to rely on 
numerical integration of Eq. H14|) . The main physics can 
nevertheless be understood in terms of the correlation- 
and retarded functions to be discussed next. 



In this section we derive and analyze the dc current 
(/) = I = Iq + Ib and the conductance G = dl/dV at fi- 
nite temperatures. We will first present the general result 
for arbitrary bias voltage, gate voltage and temperature 
to leading order in the backscattering of the two barriers 
followed by analytical approximations and a discussion of 
the results. The backscattered current is given in terms 
of the generating functional Eq. (|12(l by 



lB{x,t) = -i 



2e 



TT Sr]{x,t) 



Z'' 



(13) 



j;=0 



The actual derivation of the result is straightforward but 
lengthy. Some of the methods and intermediate results 
are presented in Appendix B. The final result for the 
current can be written as I = Iq + Ig + /^°, explicitly: 

/ = /o[l + t/'"i /"dTe*^"^^) sin[Rii(T)/2]sin(wr) 

+t/™- /dre*='i^(^)sin[Ri2(T)/2]sin(wr)], (14) 

with effective backscattering strengths [/'" = 
Em=i.2C^m' where 



U" 



{ATTtj,/h'^)e 



-Cii(O) 



E("^)'' 



A. Retarded and correlation functions 

Here, we present the results for the retarded func- 
tions and correlation functions which are carefully de- 
rived in Appendix C. In general, the retarded functions 
can be written as R^Zir) = e(r)Kfi(r) - r^^i(-r)], 

where ^'^^/(t) are the Fourier transforms of the re- 
tarded Green's functions in frequency space using a high- 
energy cut-off function exp{—\uj\/ujo). For the inter- 
acting (I) correlations at the same barriers we obtain 
^11 ("t) = ''22 (t) with 



X I e„(r) + i±^ f; 7''0a(r - 2kg) 



(15) 



Here, the interaction parameter g is introduced via 7 — 
(1 — g)/{l + g) which can be interpreted as the reflec- 
tion coefficient for an incoming charge flux traversing the 
reservoir-nanotube interface^.. For the non-local corre- 



lations we obtain r\2{T) = ''2i(''') with 



fe=0 



(2fc+l).9]. (16) 



The smeared step function is defined as 6q(t) = 
(l/7r)arctan(^) + 1/2 where a = (Ucoo)'^- The high 
energy cut-off of the theory is €q = Tiloq ~ leV which 
is the bandwidth of the SWNT. In all plots we will fix 
a = 0.001 and tip = 8 x 10^ m/s which corresponds to a 
nanotube length of L ~ 527 nm relevant for existing ex- 
periments on two-terminal ballistic transport?'^'^'^. The 
non-interacting retarded functions Rmm' (''") ^^^ obtained 
from the interacting ones by setting 5 = 1. The correla- 
tion function we decompose into a zero temperature part 
plus the finite temperature correction as 



-,I(F) 



^I(F)0 



^mm' i'^) - ^mm' i^) + ^mm' i'^) 



I(F)T^ 



(17) 



The interacting correlation functions at zero temperature 
are given as C}° (r) = C22 (t) with 



^ fc=l r=± ^ 

For the cross-terms we obtain CJ2 (t) ~ ClI (r) with 
1 --Y^ 



and the same for C^{t). In both correlation functions 
the dimensionless temperature is S = kBTtp/h. The non- 
interacting functions C^„j/ are obtained by setting g = 
1 in C]^^,. We note that all correlation- and retarded 
functions agree with the zero temperature results given 
in Rcf. 11 in the limit a — » 0. However, we note that 
a finite cut-off is crucial when doing the time integral in 
Eq. l(TH) for the case oi g = 1. 
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X E 7"= E 1^ ["' + (^ + ^(2fc + l).g)'] . 

fc=0 r=± 

(19) 

In Eqs. I|18|l and H19(l wc have dropped a r-independent 
and TTiTTi'-independent constant which does not con- 
tribute to the relevant combination Cmm'{T) — Cii(O). 
In the finite temperature part C^^t{T) the high energy 
cut-off eo can be sent to infinity (a -^ 0) as the cut-off 
is now played by the finite temperature (the result for 
finite a is presented in Appendix C). We obtain 
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fe=i 



r=± 



7r^(T + r2kg) 



sinh[7r.=,(T -|- r2kg)] 



(20) 

and the same for C2J(t). For the interference term we 
find 



Cl^2ir) 



1 
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E 

fe=0 



7- E 1- 



7rS(T + r(2fc + 1)5) 



sinh[7rS(T + r(2fe + 1)5)] 



(21) 



FIG. 3: (color online). Conductance plots: in Fig. a) we show 
the conductance as a function of bias voltage v and tempera- 
ture H both in units of the non-interacting level spacing h/t^ 
for the strongly correlated case g = 0.23. The backscatter- 
ing coefficients are taken at zero temperature as C/'" = 0.12, 
V^" = 0.1. In Fig. b) we compare the conductance at zero 
temperature for diff'erent interaction strengths: g — 0.23(red), 
g = 1 (dashed). Fig. c) is devoted to the study of the gate 
voltage dependence for g = 0.23: we have chosen [/'" — 0.12, 
17"° = O.l(red), U'" = 0.12, [/°° = -O.l(dashed), U'" = 0.12, 
17'^° =0 (points). 



B. Analytical results 

In this subsection wc provide an analytical approxima- 
tion of Ig in the regime where the bias voltage and/or 
temperature are large compared to the interacting level 
spacing h/tc where tc — Ug is the charge traversal time 
along the SWNT. In the non- interacting case g = 1, we 
can calculate Ib analytically without approximations, in- 
cluding the interference term /^°. 

Since the correlation time for the backscattering pro- 
cesses is given by fi/eV or fi/k^T the multiple reflec- 
tion terms in the retarded and correlation functions 
[Eqs. (15)-(21)] are not resolved as the traversal time tc 



is too large. We then take only the k — contribution in 
the retarded function [Eq. H15|l ] and set t to zero in the 
A; > 1 terms in the correlation functions [Eqs. p8(l and 
(|20|l ]. Note that the time t has to be considered as still 
larger than the cut-off time h/eQ. Explicitly, we consider 
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FIG. 4: Comparison of backscattered current Ib in units 
of Gofi/tpe at zero temperature (black) with the approxi- 
mate formula Eq. (I24II (dashed) showing the power-law Ib oc 
y 1-^/2 In both cases we use [/'" = U"" = 0.1 and p = 0.23. 

the regime max(ey, k^T) 3> h/tc where the incoherent 
portion P^ of the backscattered current becomes propor- 
tional to the integral 



E^ 



ex > r 

r=± 



dr 



sin(wT) 



sinh[7rS(rT -|- ia)] 



2u ■ 



(22) 



This integral can be expressed in terms of standard func- 
tions with the result (in the limit a —t 0) 



4e 



EtW) 



2 + (74^)21fcBrsinh 



/ eV 



(27r/eo)' (2nkjiT 
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\2kBT 

eV 



2nkBT 



(23) 



Note that for g = 1 the temperature dependence drops 
out and the current is only depending on the bias voltage. 
This is only true if the transmission is energy independent 
which is the case for the incoherent contribution Pg. In 
the high bias regime eV 3> k^T we obtain the power-law 
scaling on bias voltage 



2e2 



^E[W)' + ("2') 



UN21 



(27r/eo)^ 
r(2 - 7/2) 



eV 
Co 



-7/2 



-eV/eo 



(24) 



The interference contribution proportional to W^" de- 
scribes the FP-interference with two frequencies coming 
from the total charge mode with velocity Vc and three 
non-interacting modes with velocities Vp. In general, such 
integrals are not straightforward to calculate analytically 
even in the high energy regime unless we set (7=1. 

When 5 = 1, we can calculate the backscattered cur- 
rent analytically for all temperatures and bias voltages. 
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FIG. 5: Backscattered current Ib in units of Goh/Ue at 
V ~ 0.137 as a function of dimensionless temperature H. 
The continuous line is obtained from numerical integration of 
Eq. dJ, whereas the dashed line is the approximated inco- 
herent contribution oi Ib given in Eq. 12311 . As expected, they 
agree for temperatures larger than the non-interacting level 
spacing, i.e H > 1. For both curves we use [/'" = [7™ = 0.1 
and g = 0.23. 



The FP-interference contribution /^° then becomes pro- 
portional to the integral 



/p°'^~ cx ttS / dr^pf — ^-— — -, 

11 sinh(T + m^ + la)' 



(25) 



where fl — w/ttS and a ^ 0+. This integral has sim- 
ple poles (for S ^ 0) for r = — rvrS + imr — ia where 
n — 0, ±1,±2.... If fi > we can close the contour in 
the upper half of the complex plane associated with e*^"^ 
thereby picking up poles for n > 1, and in the lower 
half-plane associated with the term e~*^'^ and picking 
up poles for n < 0. All poles except the one for n — 
cancel when combining the two contributions. Adding 
the g — I contribution from Eq. l(^ (or Eq. if^ in the 
limit eg ^ oo) the result is 



/r^ 



2e2 
h 



ij 

+ Y^ lului cos(Vgi + 2Ay) 



(27rfcBT/e) . ( eVt^ 

X — r^ ; — 7^ — 7TT sm 

sinh(27rfcBTtF/n) 



(26) 



We note that temperature suppresses the FP-interference 
exponentially if kBTtY/fi ^ 1, i.e. if the temperature is 
much larger than the level spacing. 



C. Physical interpretation of dc current results 

In Section B we derived the backscattered current and 
conductance as a function of bias voltage, gate voltage 
and temperature. 

First, we discuss the incoherent contribution of Ib 
which is dominant at large bias voltages. As the en- 
ergy scale at which the system is probed exceeds the 
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interacting charge mode level spacing hvc/L, the TLL- 
correlations become apparent and our asymptotic for- 
mula Eq. (|^ applies approximately (see Figs. 4 and 5). 
At large bias voltages (and small temperatures) we ob- 
serve the characteristic power-law in Eq. H24(l . At this 
energy scale, the C/'"-term (incoherent part) does not re- 
solve the multiple reflections of the charge mode origi- 
nating from the inhomogeneity of g at the boundaries be- 
tween SWNT and reservoirs because the traversal time 
ic becomes larger than the coherence time of electron 
wave packets given by h/eV or h/k^T. In this case, only 
the k — term in retarded- and correlation functions 
contributes significantly. The strength of charge mode 
(a — 1) correlations relative to the non- interacting modes 
a = 2, 3,4 is then given as 1 — 7 where 7 = {l — g)/{l+g). 
This is apparent from the formula for the retarded func- 
tion Eq. (|15|l or in the correlation functions Eqs. H18() 
and 12UI) . This relative factor 1 — 7 is the effective TLL- 
parameter g^s at the boundary connecting a Fermi liquid 
system (metal reservoir) with an infinite TLL-systemS. 
We therefore conclude that at high voltages (or at high 
temperatures), a charge gese gets locally backscattered. 
Similar, the interference term proportional to C/'^° in- 
volves the combination (1 — 7) (1 -I- 7) = 1 — 7^ which is 
the product of two backscattering events at spatially sep- 
arated places [see Eqs. ^, ^ and ^]. The (1 + 7)- 
factor appears because the two backscattered charges can 
only interfere after the backscattered charge at the sec- 
ond barrier traverses the SWNT and has to be transmit- 
ted to the left contact with an additional factor 1 -I- 7 on 
the way. In general, the power-law behavior of transport 
can be understood as an energy-dependent renormaliza- 
tion of the bare backscattering amplitudes uj^ due to 
electron-electron interactions. It is a well known fact that 
a weak backscatterer grows strong as one approaches low 
energies, eventually going into the tunneling regime. This 
is signalled by a divergent power-law at small energies^. 
In our calculation we take into account the finite-size 
effect of the interacting region and therefore will not en- 
counter this divergence as the power-law is only valid 
above the charge mode level spacing. For sufficiently 
small bare backscattering amplitudes, the perturbative 
approach presented in this work is therefore valid at all 
energy scales. Indeed, at energies below the interacting 
level spacing, the coherence time of electron wave packets 
becomes much larger than the traversal time and even- 
tually all multiple reflections contribute. Formally this 
limit corresponds to L — > or t ^ 00 in the time inte- 
grals of Eq. H14I) where we can sum up all /c-terms and 
getting back the non-interacting functions. Therefore, 
the backscattered current is linear as V —^ 0. 

The interference contribution proportional to W^" 
shows not only a reduction of the amplitude when sweep- 
ing the bias voltage to higher values but is in principle ca- 
pable of showing FP-oscillations containing two frequen- 
cies, namely Vc/L coming from the charge mode a = 1, 
and Vp/L defining the frequency of the non- interacting 
modes a — 2,3, 4. However, the visibility of the inter- 



acting mode is in general much less pronounced than the 
non-interacting modes as can be seen in Figs. 3b) and 
3c). The reason is two-fold: First, all backscattering pro- 
cesses involve three non-interacting modes and only one 
interacting mode. Therefore, the contribution of the to- 
tal charge mode a = 1 is less pronounced. Second, the 
interacting mode contribution is further reduced by the 
smallness of gcff which enters as a prefactor in the re- 
tarded as well as correlation functions at high energies 
(i.e. small times t). At small energies (i.e. large times 
t), all multiple reflections contribute, and therefore the 
weight of charge mode oscillations increases, but then 
the charge mode behaves effectively as a non-interacting 
mode where the separation of velocities is absent. 

A finite temperature, besides triggering the TLL-effect 
[sec Eq. (|23|l ]. has an additional impact on suppressing 
the FP-oscillation amplitude which is also the case for 
the non-interacting system [see Eq. H26|) ]. This is caused 
by the temperature induced smearing of the reservoirs 
Fermi functions. But note that this suppression becomes 
exponential at large temperatures whereas the TLL-effect 
is power-law. 

In contrast to the bias voltage or temperature, the gate 
voltage does not enter as a power-law and leads essen- 
tially to a periodic modulation of If^", see Eq. (|14|l . In 
Ref. 11 it was proposed that changing the gate voltage 
Vg allows to tune the strength of ordinary FP-oscillations 
{U'^°-term) relative to the incoherent contribution (U™- 
term) which is less sensitive (through a weak gate volt- 
age dependence of g). The oscillations contained in the 
incoherent contribution is an interference effect due to a 
single impurity: the backscattered charge at the impurity 
can interfere with the momentum- conserving reflections 
due to the finite size of the interacting region. Although 
such a dependence on the gate voltage is expected, we 
note that the oscillations in the incoherent C/'"-term sur- 
vive only at small voltages eV < Uvp/Lg. This requires 
that the backscattering must be very small in order for 
the perturbative treatment to be valid. In contrast, the 
ordinary oscillations (L'^'^°-term) due to two barriers are 
more stable towards higher voltages. We will see that the 
distinction between the ordinary FP-oscillations and the 
oscillations due to the finite-size effect of the interacting 
region is much more apparent in the frequency-dependent 
shot noise. 



V. CURRENT NOISE 

The current noise S{x, lu) can be written in terms of 
the generating functional Eq. (|12|l as 



S{x,Ll!) 



'—e 

IT 



duj' 



6ri(x, uj)*5ri(x, lu') 



^r, 



(27) 



?)=0 



It is obvious from the general form of the shifted gener- 
ating functional Eq. H12|l that we can write the noise as 



backscattering is 



S{x, uj) — S'^{x, ui) + Si{x, Lu), 



(28) 



where S^{x,uj) is the noise in the absence of backscat- 
tering, and the impurity noise Si{x,uj) is the contribu- 
tion due to electron backscattering at the SWNT-metal 
reservoir interfaces. We first give the general result of 
current noise. This is followed by a discussion of the low- 
frequency noise and the Fano factor F — S/el relevant 
for existing experimentsifi. Wc then provide an analytical 
formula for the high-frequency impurity noise for general 
interaction strength g in terms of the incoherent contri- 
bution (only [/'"-contributions) which is the dominant 
source of noise at high energies. The general numerical 
evaluation including the FP-interference is presented in 
Figs. 6-9. In the non- interacting limit g = 1, we can 
calculate the noise analytically. 

The result for the current noise in the absence of 



S°{x,oj) =Gowcoth ( 



(3uj 



ReCTo(a;,2;; w), (29) 



where Re means real-part, and we have introduced the 
dimensionless conductivity^^ of the clean system without 
the backscattering cro(a;,j/;w) = {2iu!/Tr)Ri^{x,y;LLj). In 
the following, we will discuss the frequency dependence 
of the impurity noise 5*/ and refer for simplicity of the 
discussion to the situation where ^ij u^ U2 sm{VgL + 
2 Ay) — which can always be reached by tuning the 
gate voltage. We present the general result in Appendix 
B. Then, [7™ ex Y.ij ^^l^i cos(VgL -h 2A,j) will be max- 
imal for this particular gate voltage. We split the im- 
purity noise in an incoherent part plus a coherent part, 
S'/(a;,w) = Sf(x,u:) -h 5'J°(x,w). In units of GoTj/ip we 
obtain 



^/"(^'^)=-7^ E^^tl^f^^) E ko(:r,a:„,;^)P4"^(z; + r<i) 



2?iGo ^ V 2S 

r=± ^ 



2 ^ C/"j coth I — j Re o-Q (a;, a;m ; w) Im 



GQ{x,Xm\i^', 



f dT e*^"(^' sin[Rii (t)/2] (l - e'"^) cos(z;r) 



S'i°{x,uj)=-^^Re{aoix,xi;uj)*ao{x,X2;uj)} ^ coth ( ^^^^ ) ^si'^ + r'^) 



r=± 



2S 



U""" coth ( ^ ) E R-ecrn(x,x,„;w) 



X Im<^ ao{x,Xm';uj) 



/dr e^^^(^' sin[Ri2(T)/2] {dmm' - (1 - Smm')e'''^) cos{vt) 



(30) 



In Eq. (|3()|l we need the charge conductivity (Jo{x, Xm] ^) 
connecting the impurity positions Xm = ±L/2 with the 
point of measurement which we choose to be in the right 
lead, i.e. x > L/2 (the result for x in the left lead is 
easily obtained from Eq. (|31|l by xi <-^- X2 and x — > —x). 
In this case we get for the retarded function 



A. Low-frequency noise 

In this subsection we investigate the small frequency 
limit w — > of noise. 



Noise in the absence of backscattering 



BR, ^ i-^ , ^ e'^^T"^)*^ 

<(x, x„; ^) = --(!- 7) ^_^,^.,^,^ 



/'e*"(5-^)*<= j^ ^ei'^ii+^)^<'\ . (31) 



Note that i?f (x,x™;w)* = i?f (x, a;,„; -w). In Eq. ^ 
we have introduced the dimensionless frequency to = Lot-p 
and the decomposition P^ ~ Sm=i 2 ^Bm- 



We first consider Eq. (|29|l in the limit of small u where 
5°(a;,a;) ex cot'h.{f3u} / 2)uj . When U-qT ^ fiuj, we recover 
the Johnson-Nyquist noise 



5"(x,w)U^o = 2fcBrGo, 



(32) 



where we used that for small to it holds that 
coih{[3Lu/2) = 2//3a;. We also restored the units of h, 
i.e. Go = 4e^//i. In the limit of k-oT <C htu we obtain the 



quantum noise in the absence of the scatterers 

S°{x,u;)^n\uj\Go. (33) 

The full frequency dependence of cro(x, x; lu) contains in- 
terference effects of the multiple plasmon reflection inside 
the nanotube as well as interference terms depending on 
the measuring point x. Instead of elaborating the noise of 
the clean system further we concentrate on the impurity 
or shot noise, to be discussed nex1i2&. 

2. Shot noise 

At zero frequency w, the impurity noise Sj adds to the 
total noise to give 



S 



-e coth 



eV 
2fcBT 



Ib + 2kBTGB + 2kBTG, (34) 



where G = Gq + Gb with Gb = diB/dV. 
At zero temperature, the shot noise becomes S — 
-ecoih{eV/2kYiT)lB with coth(ey/2fcBT) = sgn(ey) at 
T = 0. We then finally obtain for the zero temperature 
noise at zero frequency 



el/f 



(35) 



Note that it is the electron charge e rather than the frac- 
tional charge ge in front of Ib in contrast to the infinite 
SWNT with an impurity^. 
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3. Fano factor 

Here we discuss the experimentally relevant Fano fac- 
tor F — S/el which is the ratio of the noise to the full 
shot noise el. The Fano factor is only well defined for 
the shot noise part of Eq. ^^ which is S~2kBT G. This 
Fano factor can be written in dimensionless quantities as 



coth 



eV 
2/cbT 



-'0 



2k-BT Gb 
eV Go 



1 + t 



(36) 



where we used that the total current I = Iq + Is- Note 
to be consistent with the lowest order expansion in the 
backscattering, we would have to expand the denomina- 
tor in Eq. H36|l and keep only Iq. However, this distinc- 
tion is only essential if the next order would contribute 
significantly. 



B. Analytical results of high-frequency noise 

Although an analytical solution of the time integrals in 
Eq. 1)30(1 is not possible we can estimate the general trend 
for the impurity noise at high energies. We assume that 
the temperature is close to zero, i.e. \eV\, \eV ± hLu\ ^ 
^/^cj ^B^, and obtain for the incoherent part 



St= E C/^ko(x,x,„;a;)pJ^coth 



m=l,2 



r=± 



V + ruj\\ \v + TOJ 



vr,|i-7/2 



2s ; r(2-7/2) 



2 y^ ^m coth I — j Reaoix,Xm;uj) 



(^/4) 



RCCTo 



{x,Xm;io) {sgn{Lj - v)\uj -v\^ "'^'^ + sgn{ijj + v)\u} + v\^ '^/^j 



r(2-7/2) 
lmcToix,x^;u;) cot[7r(2 - -f)/4:]i\v + Cj\^-''/^ + \v - ^\^-</^ - 2|f l^-^/^) 



(37) 



1. Analytical result for g = 1 



It is worth to examine the non-interacting limit g = 1 
in the above expression where the asymptotic approxi- 
mation Eq. (|S7|) is exact (note that the temperature de- 
pendence of the correlation functions does not contribute 
when (7 = 1, see Eq. H23|l '). We show now that we es- 
sentially get the Landauer-Biittiker result in this case, 
namely 



S" 



in, 3=1 



2h 



2_] Rn 2_\ ^'^^^ 



r=± 



eV -t- rhij 
2kBT 



(eV + rhuj) 



^-^hcOCOth (^^) h^^n^ '^"^'[^(^ + ^^)/«^] 

+ J2 ^l? cos^ [ivAx/v^] i , (38) 
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FIG. 6: The Fano factor F defined in Eq. ^ for different 
temperatures and backscattering coefficients. In a) and b) 
we show the strongly correlated case with g = 0.23 whereas 
in c) and d) we present the non- interacting case g — 1. In 
a) and c) we use C/'° = 0.12 and U"" = 0.1, in b) and d) 
we use [/'" = 0.12 and (7™ — —0.1 at the lowest tempera- 
ture. The temperatures in units of the non-interacting level 
spacing Tiv^/L are: H = O(black), H = 0.3(hght gray) and 
H = 0.7(dashed). 



where Aa; = x — L/2 > 0. In Eq. (|38|l we have written 
the noise in dimensionfuU units, and we also introduced 



reflection coefficients i?„ = Rn +Rn for barrier (1) and 
(2), respectively. They are related to C/'" by ^ J2n -^" ~ 
(7r/2)[/"^Go. This choice is motivated by the fact that 
Ib ~ — (l/2)Jo7rC/'" in the non- interacting limit (only 
incoherent contribution considered). To make contact 
with the Landauer-Biittikcr fornialism^^ we write the to- 
tal current as / = (e^//i) J2n ^"^ with T!„ = 1 — i?„ being 
the transmission coefficient for mode n. In our regime of 
small reflections (i?„ <C 1) we have / = Iq+Ib and there- 
fore 1^'^=^^ = -{e^lh)J2nRnV. The result Eq. ^ 
agrees exactly with Eq. (12) in Ref. 40 in the limit of 
zero temperature and weak backscattering. Eq. (|SS|l co- 
incides with the Landauer-Biittiker formalism only up to 
the oscillatory terms which depend on the measurement 
position X and become important if 2lo{L + Ax)/uf ^ 1. 
This oscillatory behavior for large frequencies, or equiv- 
alently, for measurement points far away from the impu- 
rities results from the beat note of finite frequency noise 
for energies which differ by ±huj. The phase difference 
acquired from the measurement point to the impurities 
and back will result in the observed interference oscilla- 
tions. 



For g = 1, we can calculate the interference con- 
tribution proportional to W^ " in cl osed form. Using 
(7r/2)C/™Go = 2{e^/h)Y.nVRiM and dimensionfuU 
units we obtain 



J 



X < Im (To (x, Xm' ; ^) cos 



2'Kk^T ^, / huj \ 

{Tiuj + reV% 



^^^"" ^ -fe sinh(2TfcBTt./n) '^"^^ [i^) ^ ^^^- ^ ER^ao(x,.„;c. 

'{nuj + reV)t 



Re ao{x, Xm' ; w) sin 



Im (To {x, Xm ; uj) cos 
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coth 
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C. Physical interpretation of shot noise results 



We first discuss the results for the low-frequency noise. 
The general result for the low-frequency noise is pre- 
sented in Eq. (|34|l . This result, valid for finite tempera- 
tures, is formally identical to that of an infinite TLI^ ex- 
cept for the important difference that a renormalization 
of the backscattered charge is absent. This fact is not at 
all trivial and has first been predicted by Ponomarenko 
and NagaosaiSS in the case of a random backscattering 



potential. Our result shows that this is also true for the 
SWNT with double barriers. One could argue that the 
high bias/temperature transport regime is sensitive to in- 
teraction (see subsection IV B) and thus a charge gese in 
front of the backscattered current appears, rather than 
the charge e. However, this conclusion is wrong. Even if 
a fractional charge is locally backscattered by the barri- 
ers, this charge cannot directly enter the leads, but will 
further get partially backscattered at the interfaces due 
to the inhomogeneity of the TLL parameter g. Summing 
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over all backscattered partial charges results in the elec- 
tron charge e. The zero-frequency noise is only sensitive 
to this integral effect as it sums up correlations over all 
times. Therefore, at low frequencies, the transport pro- 
cess is that of electrons with charge e which are backscat- 
tered by a scattering region connecting two Fermi liquid 
leads. It is interesting that this conclusion is indepen- 
dent of the bias voltage; the regime of low or high bias 
voltage is only distinguished by the power-laws of trans- 
mission. The most pronounced effect of interactions in 
low-frequency noise or Fano factor is its power-law de- 
pendence on bias voltage and temperature. 
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FIG. 8: (color online). Same parameters as in Fig. 7 but 
for a non-interacting system with g = \. In graph b) ap- 
proximately the same periodicity of FP-oscillations is found 
as in the interacting case but the frequency dependence at 
large bias voltage presented in graph c) is clearly different. 
Only much smaller ordinary FP-oscillations with frequency 
VyJL are seen. Superimposed are the additional oscillations 
depending on the measurement point. For the chosen point 
of measurement [a;| = L/2 this oscillation frequency is Vy/2L. 
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FIG. 7: (color online). Excess noise Se — S{V,u>) — S{Q,u) 
in units of Goh/tr at temperature H — 0.3 for interaction 
strength g = 0.23 using numerical integration of Eq. 1301 . The 
backscattering coefRcients are (7i" — C/2" = 0.06 and V^" — 
0.1. We have chosen the measurement point to be at one of 
the barriers. Figure a) shows the excess noise as a function of 
bias voltage eV and frequency hu) (both in units of the non- 
interacting level spacing hvp/L). In Plot b) we present the 
low frequency noise oj = as a function of bias voltage. Clear 
FP-oscillations dominated by the non-intearcting frequency 
Vf/L are seen as well as the power-law scaling with power 
1 — 7/2. In Plot c) we show the frequency-dependence of 
excess noise at large bias voltage v ~ 34. Dominant charge 
mode oscillations with frequency TiVY/2Lg are observed. 

The finite frequency impurity noise Eq. H3U|I contains 
FP-oscillations coming from all collective modes as well 
as a periodic noise suppression as a function of frequency 
with the oscillation period determined by the charge 
mode velocity Vc- At large bias voltage, this frequency 
dependence is dominated by the first term on the right- 
hand-side of Eq. (|S7|l . The periodic modulation origi- 
nates from |cro(a;,a;m; w)p where x is the point of mea- 
surement. Interestingly, this term does not depend on 
X in contrast to the terms oc coth(a)/2S) in Eq. (|37|l . 
which, however, are smaller when eV S> fiu. Therefore, 



at high bias voltage and/or low frequencies, the noise 
clearly shows the oscillations with frequency Wc/2L. This 
is indeed observed in Fig. 7. These oscillations are a 
consequence of the charge fragmentation at the SWNT- 
metal reservoir interfaces due to the inhomogeneity of g. 
These oscillations have to be distinguished from the os- 
cillations due to standard FP-interferences which contain 
the frequencies of all modes. At larger frequency uj, the 
terms oc coth((I;/2S) in Eq. (|37|) become important as 
well. They contain shot noise parts as well as thermal 
noise parts. This is clear when noting that at zero fre- 
quency these terms are the impurity dependent parts of 
2kBTGB + 2kBTG in Eq. §^. With growing frequen- 
cies UJ, these terms become sensitive to the measurement 
point X. This is clearly seen in the solution for g = 1 
presented in Eqs. (39) and (40). These oscillating fac- 
tors are a consequence of interference of electron waves 
which differ in energy by huj. On the way from the mea- 
suring point to the impurities and back these waves will 
pick up different phases which then results in the oscilla- 
tion factors, see also Refs. 39 and 23. In principle, these 
oscillations will influence the noise at every frequency, 
provided the measurement is taken far away from the 
scattering region. In reality, however, the metal contacts 
are not ballistic and therefore these results are only valid 
for a measurement point x near the barriers (within the 
inelastic mean-free path of the contacts)^. 
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FIG. 9: The excess noise Se in the same regime as in Figs. 7c) 
and 8c) but for different measurement points in the leads |x| > 
I//2. In graph a) we present the strongly interacting case 
g = 0.23 and in graph b) the non-interacting case g = 1. 
We show curves for three different measurement points: d — 
0.14(dark), d — 0.3(hght gray) and d = 0.6(dashed) where 
d is the distance from x to the nearest barrier in units of 
SWNT-length L. 



We further comment on Figs 7 and 8 which show the 
excess noise Se{V, to) = S{V, w) — S{0, uj) at low tempera- 
ture chosen to be S = 0.3, relevant for experiments. This 
definition subtracts the noise in the absence of the barri- 
ers. In the non-interacting case (Fig. 8) the most striking 
features are the clear diagonal structure in the 3D-plot 
8a) which states that the excess noise is essentially zero 
when huj > eV, see Eq. (|38|l . The small oscillations in the 
excess noise originate from the interference term Eq. H39|) 
and contain the FP-oscillations in both, bias voltage and 
frequency as well as oscillations as a function of frequency 
depending on the measurement point [see dependence on 
d in Fig. 9b)]. The multiple reflections at the bound- 
ary between nanotube and contacts are driven by the e-e 
interactions and therefore these oscillations are absent 
when g = 1. Fig. 7 shows the excess noise for a strongly 
correlated system with 5 = 0.23. The low-frequency noise 
as a function of bias voltage [Fig. 7b)] shows a power-law 
behavior with exponent 1 — 7/2 as well as some minor 
qualitative differences of the FP-interference oscillations 
compared to the non-interacting case. But the oscillation 
period is dominated by the non-interacting frequency. 
The bias window AV^ between two maxima is very well 
approximated by AV^ = h/t-pC, i.e. the voltage difference 
expected for a non-interacting system. The frequency 
dependent noise at high bias voltage is clearly different 
from the non-interacting case. Striking are the oscilla- 
tions of noise with period Aw = 7r/ic due to the charge 



flux fragmentation at the SWNT-metal reservoir inter- 
faces. They are much more pronounced than the ordinary 
FP-oscillations due to two barriers if eV ^ hv-p/L since 
P§ grows monotonically with bias voltage whereas the 
strength of /^° (showing the FP-oscillations) is bounded 
roughly by the level spacing hvp/L. At low frequencies, 
these oscillations are clearly resolved. At larger frequen- 
cies, oscillations depending on the measurement point are 
superimposed [see Fig. 9a)]. However, the charge mode 
oscillation period Tr/t^ is still very pronounced. Note, 
the excess noise is not zero anymore at large frequencies 
as this is the case in the absence of interaction. This is 
mainly due to the last term in the asymptotic formula 
Eq. H37|l which depends on the bias voltage. Therefore, 
the excess noise receives a non-zero contribution from this 
term. Note that its prefactor cot[7r(2 — 7)/4] vanishes for 
g = 1, and, as a consequence, this contribution is absent 
when g = 1. The excess noise can even get negative in 
agreement with Ref. 23. Even if S^ is not strictly zero for 
huj > eV, there is still a pronounced diagonal structure 
showing a cusp singularity when eV = fito. 

To observe the high-frequency oscillations we must at 
least be able to see the first minimum, which translates 
into Aw ~ 7r/2iF5- To have still ballistic transport we 
should be well below the mean free path of a carbon 
nanotube which, at low temperatures can exceed several 
micrometers. Using L ~ lO/im this translates into an 
estimate for the frequency of ijj/2tt ^ 100 GHz which is 
in the range of existing technology^^-^^. A relevant ex- 
tension of this setup over previously discussed systemsS^ 
is achieved through the inclusion of two impurities as 
well as through the consideration of four modes rele- 
vant for a carbon nanotube. In this system, ordinary 
FP-oscillations as well as oscillations in noise due to the 
finite length of the interacting region coexist. This al- 
lows us to extract the TLL parameter g by comparing 
the oscillations in bias voltage and as a function of fre- 
quency, i.e. building the ratio AV^/Acj ~ {2h/e)g allows 
to estimate g without referring to power-law fitting and 
without knowing of any other system parameter like the 
length L or the precise position of an impurity in the 
wireSi. We should also mention, that this ratio AV/Auj 
contains valuable information about spin-charge separa- 
tion or in general, information about different velocities 
of elementary excitations in a carbon nanotube. 



VI. CONCLUSIONS 

We have calculated and discussed in detail conduc- 
tance and finite frequency shot noise in a single-walled 
carbon nanotube (SWNT) in good contact to electron 
reservoirs using a non-equilibrium Keldysh functional in- 
tegral approach. Special focus was put on the interfer- 
ence of backscattering events off two weak impurities nat- 
urally formed at the interface between the SWNT and 
metal contacts. These so called Fabry-Perot (FP) in- 
terferences exhibit oscillations in conductance and shot 
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noise as a function of bias voltage and noise frequency 
which are dominated by the non-interacting traversal 
time tp — L/v-p rather than the interacting velocity 
tc = ^F.9, with L the SWNT-length and g the Tomonaga- 
Luttinger liquid parameter, due to two degenerate sub- 
bands in the SWNT. However, the finite frequency 
noise is in addition capable to resolve the splintering 
(momentum-conserving reflections of fractional charge) 
of the transported electrons due to the finite length of the 
interacting SWNT. This dynamics leads to oscillations in 
the frequency dependent excess noise, which, at large bias 
voltages, are dominated by a single frequency (2tc)~^, 
despite the existence of the ordinary FP-interference os- 
cillations. Therefore, shot noise measurements as a func- 
tion of bias voltage and frequency seem a decisive tool to 
distinguish the two mode velocities in the SWNT. 



APPENDIX A: CONSTRUCTION OF THE 
KELDYSH ACTION So 

The Keldysh-action So introduced in Eq. (jSl is con- 
structed in terms of equilibrium correlation functions be- 
tween the Keldysh fields. We need the following correla- 
tion function (e.g. for the 0(/>-correlation) 



'7{x,x';t) = {TKea{x,t)M^',0))o 



-({0,(x,t),<^a(x',O)})o, (Al) 



and the retarded functions 



X,x'\t) = {f KOaix,t)(t) a {x',0))o 

= -ie{t){[ea{x,t),q^aix',o)])o. (A2) 
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The expectation values are taken at equilibrium 
V — and in the absence of backscattering (ifbs = 
0). Other combinations like {TKOaix,t)(f>a{x' ,0)) = 
{TK4'a{x,t)9a{x\0)) = 0. Siucc the expectation values 
are determined by the dynamics of -ffswNT only, different 
sectors a = 1...4 do not mix, i.e. (0a(a;, i)0a'(x',O))o — 
for a y^ a' . The action Sq then can be written as 



S^ ^ \Y. I '^^ h^' I '^ [el{x,-uj),cg{x,-uj)\ Gz\x,x';cu) 

a=l 



e.{x',u) 

(j)s,{x' , to) 



(A3) 



In Eq. (|A3|I the vector 9a_{x,uj) is defined as 9a{x,uj) = [9a{x,uj),9a{x,uj)]'^ , and similar for (j)a_{x,u}), where, here, T 
means the matrix transpose. The matrix of the Green's function operator Ga(x,a;';w) — {x\Ga_{uj)\x') is constructed 
out of the equilibrium correlators and has the representation 



G'a(a;, x';uj) = 27r 



/ Cr(x,x» i?f(x,x» C'Jix,x';cj) Ri^{x,x';Lo)\ 

Rl%x',x]-uj) Rf{x',x;-uj) 

Cf{x,x'-uj) Rf{x,x'-u) Ct^{x,x'-uj) Rt*{x,x';co) 

\Rl'''{x',x;-uj) Rt'f'{x',x;-uj) / 



(A4) 



It obeys the symmetry Ga.{x,x';oj) = G'^{x' ,x;~lu) which follows from the property C^'^ {x , x' ; lu) = C^^{x' ,x;—u!), 
(and similar for 99 and i5!)(/)-correlations) which is evident from the defining Eq. IjAlfl . 

APPENDIX B: DETAILS OF NOISE CALCULATION 

In this appendix, we provide some details of the noise calculations. Starting from Eq. (|27|l we perform the functional 
derivatives and obtain to leading order in iJbs 



S{x,uj) = -^e^uj^Cr{x,x;uj) 



Q r 



rr'^zk 



J2 rr'Ux,u) [f:,{x,^Lo)~e-^-'rr'A^,~u) 
I 



.»[^™„(t)-^"'„(0)] 



, (Bl) 



where f^{x,uj) — uj[Ci^{xm,x;u}) ± |i?J*(x,Xm; — oj)] and the phase operators are A^^,g{t) — 9^^ + s92 



Further, we intro- 



and ^^„, = [V^L + 2A„)(m - m') 
duced the index n = ijs and the abbreviation (...)o = 
fl^ J'I?[6'^0^]... exp[iS'o]. We assume now that the cur- 
rent is measured in the right lead (the result for x in the 
left lead is easily obtained by xi <-^- X2 and x — > —x), 
i.e. X > L/2. The general expression for the retarded 
function in that case reads (see Appendix C) 



i?l {x,Xm;Uj) 



(1-7) 



,(l-SgL)t^ 



- 76 



Ai+^)t. 



(B2) 



Note that Ry{x,x. 



— Ri^{x,Xm',—^)- The cor- 
relation functions are related to the retarded functions 
via the fluctuation dissipation theorem Cf^{x,Xm;u!) = 
Cf^{xm,x;uj) = — coth{f3uj/2)ImRf^{x,Xm',^), where 
Im denotes imaginary part. The expectation value in 
Eq. IJB1|) is of the general form (e, e' = +7 — ) 

(B3) 

where we used that the action 5*0 is quadratic in the 
bosonic fields 9^ and cj)^ which allows to perform the 
average in the exponent. The correlator in the exponent 
of the right-hand-side (RHS) of Eq. (|B3|I always (i.e. for 
general n) contains a sum of three non-interacting modes 
a = 2, 3, 4 due to spin and subband degeneracy plus one 
interacting mode of the total charge a = 1: 



(A;;,„(i')±<'nr) 



^ [0:^(1') ± 9:^,(1"))'+ 6,, Y. (^a™(i') ± oL'{t") 



a— 3 



-(l-'5.,)E(<^-n(*')±0a™'(O) ) ■ (B4) 

/ 



a— 3 



15 



ings which, however, leads to the same correlation func- 
tions and only the scattering phases hidden in V^" dis- 
tinguish the different processes. Also note that the cor- 
relation functions do not depend on the spin direction s. 
We find [e.g or the ^-fields (similar for 0-fields)] 



Note that the correlator above depends on the different 
processes of interband i 7^ j or intraband i — j scatter- 



J 



(0^„(i)±C'(O) 



±2C! 



'W + 2Cff™™(0) 



±^ir'-r) KUn.'it)-KU'rai-t) 



±-{r + r') 



R. 



.'it)+KU'^n{-t) , (B5) 



where we have simplified the notation for clarity of 
the presentation: C^^„,rn'i^) - C^^{xm,Xm';t) and 
Rl1^^,it) = Rl%x„,,x„,,-t). We note that only the 
- sign in Eq. IJB5(I contributes as exp[— (l/2)((0JJ,„(i) + 
^am'(O))^)] « exp(-oo) due to the first line of the RHS 
in Eq. (jBSp . This is a direct consequence of particle 



conservation since the -I- sign option comes from terms 
like {ip^ip jiip^ip r) Q . The general form of the noise de- 
scribed by Eq. (|Bip can be split into a sum of a clean 
limit with no backscattering corrections S''\x, w) and the 
backscattered correction Si(x, uj) which we refer to as the 
impurity noise. The clean limit can be written in a more 
standard formS^ using the relation between retarded and 
correlation function 



S'^{x,uj) — Go^coth 



f3uj 



'ReaQ{x,x;uj), (B6) 



where Re means real-part and we have introduced the 
dimensionless conductivit}i2i of the clean system with- 
out the backscattering ao{x,y;uj) = {2iuj / tt) Rf^ {x , y; uj) . 
The impurity contribution to the noise Si{x,lu) can be 
calculated using Eq. (|BH) and Eqs. l|54|) and (|55)l . We 
obtain after some calculation the general result valid for 
all temperatures, frequencies, gate voltages and to lead- 
ing order in the backscattering Hamiltonian Hhs (in units 
of Goh/U) 



Si{x,uj) = -2coth I — j Yui^Reao{x,Xm;^) 



X Im 



ao{x, x„,;uj) f dr e*^"(^) sin[Rii(T)/2] (l - e*"^) cos(wr) 

^^gCii(r) sin[Rii(r)/2] sin[(w + rC:i)T] 



"2-^ '"l^"^^'^™''^-" Z^coth 

m r— ± 



2S 
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- coth f — j ^ ReCTo(x,Xm;w) 

mm' 

X lni{aa{x,x^,-uj) f dr e^^'^^^ sin[Ri2(T)/2] {d„un' - (1 - S„un')e'^^) [C/'° cos(z;r) - 1^™(1 - 2Si^,) sm{vT)] 

+ ^Relaoix,xi;Luy(Toix,X2;uj) /"dre'^^^^^) cos[Ri2(|r|)/2] e*"^[t7'=°cos(wT) - T/™sin(wT)] | . (B7) 



In Eq. (im we have used Timm'ir) = ^Lm'('^) + 
3i?J^j„/(T) and a similar definition holds for Cmm'{T). 
The superscripts I and F denote interacting and 
free, respectively. The interacting functions are 9i9i- 
correlations whereas the free functions come from cor- 
relations of the non-interacting modes a = 2,3,4. In the 
main text we give the slightly more compact result for 
the case where the coherent (FP)-contribution is maxi- 
mum, i.e. when \U'^°\ is maximum as a function of gate 
voltage. In Eq. IjBTII we have introduced V^" which is just 
C/™ with cos(FgL + 2Ay) replaced by sin(VgL + 2Ay). 
We see that at finite frequency w, the impurity noise be- 
comes sensitive to the real and imaginary parts of the 
conductance (Jo{x, Xm] ^) which contains the multiple re- 



J 



r 



flections of the charge mode a = 1 at the inhomogeneity 
of g where the S WNT is connected to the non- interacting 
reservoirs. The complete noise as a function of frequency 
and bias voltage is therefore a complicated superposition 
of ordinary FP-oscillations described by the time inte- 
grals which are influenced by both, voltage and frequency 
and exhibit by all four modes whereas the additional fre- 
quency response due to UQ{x,Xm\'^) is only sensitive to 
the total charge mode a = 1. For clarity, we give here 
the explicit form of real and imaginary-parts of the re- 
tarded Green's function connecting the two barriers with 
the measurement point x (assumed to be in the right 
lead). For the retarded function with .t,„ = —L/2 we 
obtain from Eq. HB2p 



i?i (a;,xi;a;) 



1-7' 



2w 1 - 272 cos(2cjic) + 7' 



^ {sin [ujU{g + d)] + 7^ sin [ujU{g - d)] ] 



— I ■ 



1-7' 



2a;l-272cos(2wic)+7'' 
For the retarded function with x„i — +L/2 we obtain 

1-7 



{cos [ujU{g + d)] - 7^ cos [ujU{g - d)]] . (B8) 



Rl%X,X2]L0) 



2oj 1 - 272 cos(2wic) + 7^ 



{sin [iot-pd) [1 + 7(1 — 7) cos(2cL'ic) — 7'^] +7(1 + 7) cos [uut-pd) sin(2wtc)} 



1-7 



2uj1~ 272 cos(2wtc) + 7' 



{cos {ijjUd) [1 + 7(1 - 7) cos{2uotc) - 7^] - 7(1 + 7) sin [uoUd) sin(2wtc)} • (B9) 



In Eqs. IIB8(l and (|B9(I we have introduced the distance 
from the measurement point x to the nearest metal 
contact-SWNT interface (here X2 = L/2) in units of the 
length L of the nanotube, i.e. d ^ {x — X2)/ L. 



For completeness, we also show a more direct way to 
obtain the low-frequency noise Eq. H34|l starting from 
Eq. (JB1(I using a low- frequency expansion in uj. For 
this we consider the term J2rr' ''"''' fm{^T^)[fm{^i^^) ^ 
exp(— jtt't)/^,(a;, —uj)] in the limit a; — *■ 0. To proceed in 
the evaluation we flrst make a straightforward expansion 
of the above expression. We write the function f^{x,uj) 
in terms of real- and imaginary parts of the retarded 



I 

function Rf^{x,Xm',i^) as 



^-coth^^) ) Imi?f(x,x™;w) 



+ -rReRi {x,x,n;uj) 



(BIO) 



We now use the low-frequency behavior of the real- and 
imaginary part of the retarded functions in Eq. (jB10|) 
which we express as 

Rei?f (x,a;,„;u;) ^ Rl°{x,x,n) + 0{uj^) (Bll) 

and 



Imi?i {x,Xm]Uj) 



2lo 



+ R^^Hx,x^P, (B12) 
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where R^^{x,Xm) and i Rf^{x,Xm) are the zeroth-order 
and Ist-order expansion coefficient of Rf^{x,Xjn',u!), re- 
spectively. They depend on x, Xm and g but we find that 
these terms will not contribute to the zero- frequency limit 
of noise. As a consequence, the low-frequency noise is in- 
dependent on the position of measurement. In this limit, 
we obtain for the frequency dependent part of Eq. (|B1|) 



rr' f^{x, uj) f^{x, -ui) - exp{-iujt)f^, {x, -w) 



-^(.fcBT)^t^-(l --')(£' 

-f 27r [k^Tf rr' {R{Hx, x^) - i?f (x, x^')) 

-i-k^T {r'Rl^[x,x„i) - rRl°{x,x„r>)) 



-0{uj). 



(B13) 



By performing the sum over rr' as well as the sum over 
mm' we find that the time integral in Eq. (|Bip yields 
zero for the terms associated with the 1/lu contributions 
in Eq. IJB13|) . Therefore, the limit w — > is well de- 
fined. Note also that the voltage term cos{eV t+ip'^^,) = 
cos{eVt) cos{ip'^j^,) — sin{eVt)sin{(p^^,) where, due to 
the symmetry in the sum over mm' , only the cosine-terms 
contribute. Using Eq. 1|B13|I in Eq. (|Bip for the noise we 
obtain the low- frequency noise presented in Eq. H34|l . 



APPENDIX C: DERIVATION OF RETARDED 
AND CORRELATION FUNCTIONS 

In this appendix, we outline the derivation of the re- 
tarded Green's functions which are calculated in the equi- 
librium system and without the backscattering (_ffbs=0). 
We choose to calculate the temperature Green's function 
first and then rotate back to real time (Wick rotation) 
which gives us the retarded function. 

We start with deriving the action from the Hamiltonian 

-fJSWNT 

>C(n„0J= fdxJ2 [^a{x)ea{x)-H,^,^riIla,ea] 
'' a 

(CI) 
The action is defined as the time integrated Lagrangian 
S — J dt£{t). We then change to imaginary time r — 
it and introduce the Euclidean action by the standard 
identification ~Se — iS{it — > r). This immediately gives 



«-4i: 



/3 



dr / dx 



2 V 



9a 



:{d^eaf 



(C2) 



To calculate time-ordered correlation functions 
(TA(x)i3(x')) where A and B are any function of 
operators Oa and 0^ we can use the functional integral 
approach 



SB(e„,0„) 



(fi(x)i3(x')) - |n fv[q^aOa]A{^)B{^')e' 

a •' 

(C3) 
where x = {x, r) and Z = ll^J V[(t)aOa\ e-^^^^-"^") is the 
partition function. Here, since the bosonic fields are her- 
mitesch, the functional integral is over real- valued fields. 
If the operators A and B are only functions of one of 
the field- type, i.e. only a function of either 6a or (pa, we 
can integrate out the other variable to get an effective 
action which only depends on one of the variables. To 
do this we use the result for Gaussian integration over 
multidimensional real variables Xi, i = 1, ..., N 



n 



dxi 



5 J^ XiAijXj+Y^ \iXi 



/ e *^ 
VdetA 



(C4) 
where det means the determinant of the real symmet- 
ric and positive definite matrix Aij. We first derive 
the action for the 6a— fields. Using partial integra- 
tion in the action Se we can bring the functional in- 
tegral Eq. (|C3|) to the form of Eq. IIG4|l with identify- 



ing Xi = {i/TT)dxdr6a and the matrix elements Aij = 
— {vF/7T)d'^S{x — x')6{t — t'). Note that the determinant 
will cancel with the similar integration procedures in the 
partition function. Therefore, for calculating correlation 
functions the explicit calculation of the determinant in 
Eq. IIG4p is not needed. We then obtain the effective 
action for the 6„— fields 






1 

in 



E 



dr i dx 



-{dx6af 
9a 



1 



9aVa 



idr6af 



(C5) 
Similarly, for the effective action of 0a— fields we obtain 






dr / dx 



Va9a{dAaf + —{d^<Paf 
Va 



(C6) 
We can now use the effective actions to calculate the 
correlation functions (r0a(x)6'a(x')) and (T(/)a(x)0a(x')), 
respectively. 

The relation between functional integrals and time- 
ordered correlation functions Eq. I|('3|l (e.g. for 0a-fields) 



(re,(x)0a(x')) = ^ / I?0a0a(x)0,(x')e-^-(''" 



can also be written as 

(m(x)0a(x')) = .,. :.,. .. z', 

where 



A=0 



5A(x)(5A(x') 



(C7) 



(C8) 



(C9) 
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is the generating functional for the ^a-field. Writ- 
ing the action S£ as a biUnear form 5*^° = 
(l/2)/dx/dx'6'a(x)[G^^]-i(x,x')6'a(x') and using the 
result for Gaussian integration Eq. (|C4p we conclude that 
{T9a{x)9a{x')) — (^^^(x, x') accompanied with the oper- 
ator statement 



G 



Gf = 1. 



(CIO) 



Using that the inverse Green's function operator G ^ is 
local in (imaginary)time and space, i.e. (x|G^^|x') — 
D{x.)S{x — x'), leads us to the differential equation for 
the Green's function in imaginary time 

D{x, t)G{x, t; x',t') = 6{x - x')6{t ~ t'). (Cll) 

Eq. IjClljl clearly shows that the Green's function is sym- 
metric in X and x'. Explicitly, we obtain the differential- 
operators £)^^x) = -l.[-l-d^, + d^^d,] and L)f'^(x) = 

— -[^d^ + dxVagadx]- Due to the inhomogeneity of the 
charge mode a = 1, its Green's functions are not trans- 
lational invariant. The time translation invariance al- 
though still holds, i.e. G(x, r; x' ,t') = G{x, x'; t—t') and 
we can transform to frequency space using the Fourier 
expansion for boson Matsubara Green's functions 



G{x,x';t) = ^^^e' 



'G{x,x';ujn), 



(C12) 



where the sum is over the Matsubara frequencies w„ = 
27Tn/(3 with n = 0, ±1, ±2. We then obtain a differential 
equation for the Fourier component of the Matsubara 
Green's function 

-^ivl-d^^dx) Gi'{x,x';iJn)=TTSix-x'), (C13) 

9aVa 9a J 

and 

(^u;l-dxVagad.^ Gf{x,x';u;n) = ttS{x-x'). (C14) 

The solution to these partial differential equations for 
the case \x'\ < L/2 can be found by the ansatz^^(for the 
6'0-contribution) 

Gr(x,a;';c^„) = 



Sel'^-k/fa ^ c'g-h^k^a; _L/2 < x < x' < L / 2 
D^\^r.\xlv^ + £:g-|^„b/t,„. _^/2 <x' <x<Ll2 

P^-\u^r,\x/v,. X>L/2 

(C15) 

These solutions satisfy the boundary condition 
Gl%±oo,x';uJn) = 0. The coefficients A ~ F are 
functions of x' and ujn and can be found from the 
following boundary conditions 

• Gff{x,x';ujn) is continuous everywhere. 

• ^^^^'^xGff{x,x';uJn) is continuous at x = L/2, —L/2. 

• — ^^j^dxGff{x, x'; LOn) has a step of height tt at x = a 



These three conditions lead to the following set of equa- 
tions which can be used to determine all constants A — F: 

F = _C)el'^"l-^(i+9")/^"F -I- iJe'""!"^'"'""^"'/^"'' 
G + Be^l'^^l^'A-o ^ E + £)e2|'^"l^'/"a 

g^A = Bel""l^(^"9a)/2«F _ (^gh„|L(l+g„)/2i>P 

q^P = £;gl'^i.|i(l-9a)/2DF _ £)g|w„|L(l+5a)/2t)F 

B-L»+(S-G)e~2|a;„|:E7t,, 



QaT^ „-\lj„\x' /v„ 



The retarded Green's function R^^{x,t;x',t') = 
-iQ{t~t'){[ea.{x, t), Oaix', t')]) is obtained from the Mat- 
subara Green's function via the analytic continuation 



Rl\x,x';iu)^-Gl'{x,x';u;r.)\ 



(C16) 



with 6 = 0^. The analytic continuation is performed 
from the positiv imaginary axis where the function is 
G^^ {x , x' ; oJn) with tJn > to just above the real axis 
where it equals R^^ {x , x' ; uj) . This amounts to the re- 
placement \u!n\ -^ — iw -(- (5 in G^^{x, a;'; w„) to obtain the 
retarded function R^^{x,x';uj). The function for the (f>a- 
fields can be obtained from the solution for the 6'a-fields 
by the substitution ga -^ ^/ga which is evident from the 
differential equations Eqs. IIC13I) and (|C14|I . For the 
retarded functions i?^f,„(cj) and R^f,Juj) we need the 

alii*- ^ a|ll^ ' 

solutions in Eq. (|C15p with —L/2 < x < L/2. Since the 
Green's function is continuous everywhere we get also the 
correct solution at the boundary where x' — zLL/2. We 
obtain for general |a;,a;'| < L/2 and for the interacting 
mode a = 1 



Ri [X, X ;u!) 



''''"'9 (iat^\x~x'\/L I 



7 



2lu 



p — i2ujtc ^2 



X y^ \^-iuit^[l-r(x+x') / L] _|_ ^iu)t^r(x-x')/L 



, (C17) 



where Co — uj + i5, and tc — Lg/vp. We further introduced 
7 = (1 — g)/{l + g) which can be interpreted as the re- 
flection coefficient for an incoming current flux traversing 
the reservoir-nanotube interfacc^^(i.e. the inhomogene- 
ity of 5). We also need the retarded functions in real time 
which we get by Fourier transforming Eq. ljC17p . Using 
(1 — -y2gi2wtc-j-i _ ^00^^ ^2fc gi2fci^t<, g^j^j ^ high-energy 

cut-off function e^''^'/'^^ we obtain for x = x' = ±L/2 



4Ut) = -T;il-l) 



^{&^oit) + 



1 + 7 

7 



^7''e^„(i-2A;t,) , (C18) 



k=l 



and for the cross-terms x — L/2{—L/2), x' — —L/2{L/2) 
describing the FP-interference we obtain 

00 
<i2W = -f(l-7')E^''©-[t-(2fc + l)ic]. (C19) 

A:=0 
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The smeared step-function is defined as Qujoit) = 
(I/tt) arctan(aJoi) + 1/2. If we keep the cut-off finite 
the correct retarded function is obtained by the com- 
bination i?ff ,(i) =6'(i)[r^f 

l\mm' y ' ^ / L 1 



Note 



,(t)-rf, , (-t)\. 
that the retarded Green's functions are temperature in- 
dependent. The temperature dependence is completely 
contained in the correlation functions to be derived next. 
First note that we are dealing with equilibrium proper- 
ties, and therefore the correlation function is connected 
to the retarded function via the fluctuation-dissipation 
theorem 



O. 



a\r 



(C20) 
We give here the results for the 06'-correlations. The cor- 
responding results for the (/)(/)-correlations are obtained 
by the replacement: ga -^ ^/9a- We split the tempera- 
ture dependence in a T = part plus the temperature 
corrections, Cff„„,(t) = C'^ 
that we can decompose 



'it) + C:Zn'it)- Note 



coth ^— = IH 5-. 

2 / l~e-0^ 



For positive frequencies lu we can write 1/(1 



(C21) 



-I3uj 



) as 

a geometric series valid for all temperatures. We then 
obtain coth(/3cj/2) = 1 -f 2X;,T=o ^"''"^"^^^ ^^ich cor- 
responds to the two terms contributing either at zero 
temperature coth(/3ct;/2) — 1 or to the finite temper- 
ature corrections coth(/3w/2) = 2X;r=o '^'^"^"'''^^- ^^ 
zero temperature we then obtain 



C^aTl2W 



](l-7^) 



J2 7'' E l'^ K' + (* + '^(^fc + l)tcf] , (C22) 



fc=0 



■r=± 



and 



Note in Eqs. IJC22|) and IIC23p we omitted a time and 
space-independent constant which does not contribute to 
the relevant combination Cf^„„,(i) - Cf^„„^(0). For the 

finite temperature correction C^f^ we obtain 



cZli*) 



X Vt^'-^Viu 



fe=0 



r=± 



1 

LUoP 



■ t+rtc{2k+l] 
' 



r 1 



1 



(C24) 



In the finite temperature correlation functions it is al- 
lowed to perform the limit loq ^ 00 as the finite temper- 
ature plays the role of the cut-off. This is true as long as 
k^T is small compared to the high-energy cut-off eq. Af- 
ter doing so, we can use that |r(l -I- ix)\'^ = ttx/ sinh(7ra;) 
for X real to obtain the simpler form 



c:tm - 



1-7' 



00 

x5:,-Ein 

fc=0 r=± 



TT{t + r{2k + l)tc) 
/3sinh[7r(t + r(2fc + l)tc)//3] 



(C25) 



For the autocorrelation functions {m = ni') we obtain 



L-^l«,,ri(tj 



1-7 



2 0° 



In 



Tit 



1-7 ^ 
47 



E^^'Ei- 



/3sinh(7ri//3)_ 

7r(i -I- r2ktc) 



fc=i 



r=± 



/3sinh[7r(i + r2fctc)//3] 



(C26) 



^.TnW--V^"(""+^') 



1 + 7 
27 



E '^"' E 1^ K"' + (* + r2kt,f] I . (C23) 



k=l r=± 



The frequency representation of the retarded function 
Ri^{x,Xm',(^) given in Eq. (jB2|) where the measurement 
point X > L/2 is chosen to be in the right lead and 
Xm — ±L/2 can also be obtained from Eq. I|C15|) in the 
regime x > L/2. 
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